Elastically-Driven Linker Aggregation between Two Semi-Flexible Polyelectrolytes 
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The behavior of mobile linkers connecting two semi-flexible charged polymers, such as polyvalent 
counterions connecting DNA or F-actin chains, is studied theoretically. The chain bending rigidity 
induces an effective repulsion between linkers at large distances while the inter-chain electrostatic 
repulsion leads to an effective short range inter-linker attraction. We find a rounded phase transition 
from a dilute linker gas where the chains form large loops between linkers to a dense disordered linker 
fluid connecting parallel chains. The onset of chain pairing occurs within the rounded transition. 



Highly charged chains of the same sign can be linked 
together by mobile polyvalent ions of the opposite sign 
(polyvalent counterions) . This phenomenon has been ob- 
served for a wide variety of systems, including solutions 
of DNA HI, F-actin |-§ and polystyrene sulfonate §. 
For flexible chains, polyvalent counterions can induce col- 
lapse of the chains into a globular compact structure ; 
alternatively, one can argue that the flexible chains medi- 
ate attractions between polyvalent counterions that cause 
them to aggregate. For rigid rods, on the other hand, the 
polyvalent counterions always repel each other along the 
axis of the rods, and the attraction between rods is at- 
tributed to ion-ion correlations among rods |^|,^) . 
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FIG. 1. Schematic view of two chains connected by linkers at 
(a) low linker densities (the dilute regime); (b) high densities (the 
disordered "railway track" structure) and (c) intermediate densi- 
ties (the crossover regime). One chain is in the y = plane while 
the other is in the y = d plane. 



In this paper, we will study the intermediate case of 
semi-flexible chains, appropriate to biologically impor- 
tant polymers such as DNA. It has been shown theoret- 
ically that counterion correlations can modify the bend- 
ing rigidity of a single semi- flexible chain |l0| |l2| and can 



even render the chain unstable to collapse jLl|. Here, we 
consider two semiflexible chains from a different point of 
view: instead of studying how counterions modify the 
effective interactions between monomers on chains, we 
examine how chain flexibility modifies the effective inter- 
action between generalized linkers, which could be simple 
polyvalent counterions H or weakly-binding (crosslink- 
ing or bundling) proteins Q. Alternatively, the link- 
ers could represent hydrogen bonds connecting the two 
strands of a DNA double helix undergoing denaturation 
jT3yri] |, We use this effective interaction to study the 
many-body statistical mechanics of linkers. A similar 
approach has been fruitful for understanding behavior of 
proteins that link together elastic membranes [ p>[]l6| ] . 

Our calculations yield three main results. First, we 
find that the chain-mediated interactions between link- 
ers are non-monotonic. At large linker separations the 
chain bending elasticity leads to a long-ranged repulsion, 
while at short distances the electrostatic repulsion be- 
tween chains leads to a short-ranged attraction between 
linkers. Consequently, there is a repulsive barrier in the 
interaction between two linkers at intermediate separa- 
tions. Second, the unusual shape of the effective poten- 
tial leads to interesting phase behavior in the many- linker 
system. Since the two-chain system is one-dimensional, 
there is no true phase transition instead, we find a 
rounded transition from a dilute phase of linkers where 
the chains form large loops to a dense one where the 
chains are parallel and close together. The rounded tran- 
sition is accompanied by large fluctuations in the spacing 
of linkers. This is reminiscent of large fluctuations in the 
separation of membranes near their binding/unbinding 
transition ]is| . Third, we find that the onset of aggrega- 
tion of single chains into pairs occurs within this rounded 
transition. Our results suggest that in many-chain sys- 
tems, the rounded transition might be found near the 
onset of bundling. 

Our model consists of two charged semi-flexible chains 
held together by a series of linkers (Fig. |l|). Each chain 
carries a negative linear charge density A = —e^/lg where 
e is the electronic charge, Ib = e 2 /eksT is the Bjerrum 
length, e is the solution dielectric constant and feT is 
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the thermal energy. Here, A and £ are the effective lin- 
ear charge density and the corresponding, dimensionless, 
Manning- Oosawa parameter HJljJ, respectively. In the 
presence of salt at concentration q,, the Debye-Hiickel 
screening length = I/^/EttTbCI) characterizes the ex- 
ponential decay of electrostatic interactions. We assume 
that a linker of size b constrains the two chains of radius 
r s at a fixed separation d = 2r s + b. 

The first step is to calculate the effective interaction 
between two linkers separated by a distance I. This en- 
ergy can be estimated by considering the total energy 
of a periodic array of linkers with separation distance I 
along the x-axis. For simplicity, we assume that each 
chain is restricted to a plane parallel to the xz-plane. 
One chain is located in the y = plane and follows the 
curve {x,y, z} = {x,0, z(x)} where z(±l/2) — and 
z(Q) = h/2. The second chain is located in the y = d 
plane and follows the curve {x,y, z} = {x, d, — z(x)}. A 
periodic trial function is assumed for the two chains 



z(x) = — cos (ttx/1) 
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where h is related to the crosslink angle 9 through 
tan(6>/2) = irh/21 and will be determined variationally. 
It should be emphasized that although the exact shape 
of the trial function affects the numerical values it should 
not have a qualitative effect on the results. 
The electrostatic energy per linker is given by 
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where rf 2 = (x — x') 2 + d 2 + (z(x) + z(x')) 2 and z x = 
dz/dx. In Eq. ^|and in the following discussion all ener- 
gies are expressed in terms of the thermal energy fe^T. 
Note that if two straight chains cross at an angle 9, the 
electrostatic energy is [fill 
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favoring a crossing angle of 9 = 90° and providing a basic 
energy scale in the problem. 

This interaction competes with the bending energy EG] 
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which favors small crossing angles. The chain persistence 
length l p characterizes the decay of correlations in the 
chain orientation. Any dependence of bending rigidity on 
the composition of the surrounding solution |n^-[l2|,^l| is 
implicit in l p . 

The total free energy per linker f(l;h) = E c \ + -E^end 
can be calculated numerically and minimized with re- 
spect to h p2]. Typical interaction energies A/(Z) = 



/(0 — as we ^ as the behavior of the crossing angle 

9, are depicted in Fig. |[ Note that f(l — > oo) = T. 




FIG. 2. (a) Interaction free energy per crosslink A/(i) = f(l)—T 
and (b) crossing angle 8 as functions of the reduced distance kI be- 
tween linkers. Note that A/(/) changes sign near lo (see text), 
and that 9 becomes nonzero at Z* (denoted by an arrow for the 
solid curve). The curves correspond to different values of the chain 
persistence length l p = 10 A (dots), l p = 50 A (long dashes) and 



l p = 200A (solid curve). We use d = 25A, £ : 
l B = 7k. 
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The competition between inter-chain electrostatic re- 
pulsion and intra-chain bending rigidity leads to an effec- 
tive attraction between linkers at short separations and 
an effective repulsion between linkers at large separa- 
tions. At short separations, the bending energy domi- 
nates over the electrostatic energy and the chains remain 
parallel to each other (crossing angle 9 = 0). Since the 
chains are straight, the bending energy does not other- 
wise affect the linker interaction. Furthermore, the total 
electrostatic repulsion between chains is fixed and does 
not depend on the position of the linkers. As the den- 
sity of linkers increases, the energy per linker decreases. 
The effective interaction can be estimated from the elec- 
trostatic energy of two parallel chains (T^|. It is at- 
tractive and depends linearly on the linker separation: 
A/ s hort — — r(7 — l)/h where nl Q = ir/e Kd K (Kd) and 
Kq is the zero-order modified Bessel function. The origin 
of the attraction can also be understood by considering 
only two linkers: by joining to form a single junction 
instead of two separate junctions, the electrostatic repul- 
sion between the chains is reduced by roughly the junc- 
tion energy T. Thus, the inter-chain repulsion leads to 
an effective inter-linker attraction. 

At very large inter-linker separations, the electrostatic 
repulsion dominates the bending energy and the cross- 
ing angles saturate to 9 — 90°. Since the chains are 
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perpendicular at the junctions, the electrostatic repul- 
sion does not otherwise affect the interaction between 
linkers. Compressing the linkers together costs bend- 
ing energy, leading to a long-range repulsion of the form 
A/i on g ~ alp /I, where a is a constant of order unity. We 
estimate a ~ tt/^/2 by assuming constant curvature of 
the two chains in between junctions and minimizing with 
respect to 9. The crossover between the two regimes oc- 
curs at a separation I* where |A/i ong | — |A/ sri0 rt|- For 
lp ~ r^o (as is the case for the parameters of Fig. |^) 
we obtain I* ~ Iq + al p /T. For stiffer chains, we find 
I* w y/al p l /T. 



Zi = 
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The final inter-linker potential we use is 
v(l) = Af(l) + cki(l/l p ) 



(5) 



where the second term estimates the entropy loss of two 
random walks of step size l p constrained to cross after 
the same number of steps at a distance I from each other 
jL3j. We use c — 3/2 which corresponds to ideal chains 
in three dimensions. A recent estimate which includes 
excluded volume interactions between different parts of 
the chains yields c ~ 2.1 — 2.2 0|. This term is relevant 
only at large distances I 3> l p . In addition, the linker size 
b provides a lower cut-off on the inter-linker separation I. 
This corresponds to a hard-core repulsion or, in the case 
of DNA denaturation, to the size of one base pair. Other 
direct interactions between linkers such as the Coulomb 
repulsion between polyvalent counterions can also be in- 
cluded p2] but do not affect our main results. 

With the effective linker interaction now in hand, we 
can study a one-dimensional fluid of N interacting link- 
ers located along a straight line of length L at < x\ < 
X2 < ■ ■ ■ < xm < L. At low densities (large separations 
between linkers) one might expect a crystal of linkers 
even in this one-dimensional system because the effec- 
tive repulsion decays as l/l. However, in this regime, the 
crossing angle is 9 = 90° and the chain configuration on 
one side of the crossing junction is decoupled from the 
chain configuration on the other side. As a result, the 
linker interaction is only a nearest neighbor interaction; 
that is, the interaction decays as l/l, but only as far as 
the nearest linker. 

For nearest neighbor interactions, one can calculate the 
partition function directly from the inter-linker potential 
v(lj = Xj+i—Xj). Note that at higher densities, the cross- 
ing angle becomes less than 9 = 90° and further-neighbor 
interactions develop. Nevertheless, at those shorter sepa- 
rations, the interaction no longer has the form l/l, so the 
further-neighbor interactions cannot lead to new phase 
transitions, although they might affect quantitative re- 
sults. We therefore neglect further-neighbor interactions 
at all densities. 

It is convenient to use the Gibbs ensemble in which 
the natural variables are the number of linkers N and 
the one-dimensional pressure P E^j. In the thermody- 
namic limit the Gibbs free energy is G = —NhiZi where 



and At is the thermal wave length. The chemical po- 
tential is p = —hxZi and the linker density p = N/L 
can be calculated from the average linker separation 
1/p = I = — d(ln Z\)/dP. The fluctuations in the sepa- 
rations Al 2 = (l — V) 2 can be similarly extracted. 




FIG. 3. One-dimensional linker density p as a function of the 
chemical potential /i and (inset) the relative fluctuation in the in- 
ter-linker separation Al/l as a function of p. The parameters are 
the same as for the solid curve in Fig. ^| for which T ~ 11.8; also, 
b = lA and X T = lA. 

The typical behavior we obtain is shown in Fig. ||, 
where we plot the linker density p as a function of chemi- 
cal potential p (shifted by T — 2e where £o is the adsorp- 
tion energy of a linker onto a chain) . Note that adsorp- 
tion energies and a constant T term are not included in 
the definition of p. At low p the linker density is low and 
the chains form large loops between the linkers (Fig. |l]a) . 
At high p, on the other hand, the system resembles a 
disordered railway track where the chains are rails and 
linkers are ties (Fig. |l]b). 

Between these two limits, we find that at a certain 
value of the chemical potential p* < — T the linkers ag- 
gregate together and the density jumps sharply. The 
slope at the jump appears vertical because dp/ dp is large. 
Within our model, there is no true phase transition [JL7J. 
However, this jump in p can be viewed as a rounded tran- 
sition from a dilute to a dense phase of linkers. In the 
vicinity of the rounded transition there are large fluc- 
tuations in the inter-linker separation, as shown in the 
inset to Fig. || These fluctuations correspond to clus- 
ters of linkers that increase in size as the system crosses 
through the rounded transition. The crossover regime is 
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limited to a narrow region around p* but spreads over a 
considerable range of densities (0 < p < 0.2 in the fig- 
ure). Finally, we find that as the density increases near 
the rounded transition, the short-ranged attraction be- 
tween linkers sets in and the free energy drops rapidly. 
As a result, paired chains become favorable compared to 
isolated chains. Thus, the transition from single chains 
to pairs occurs near the rounded transition (/x ~ /i*). 
The density p* at which chain pairing occurs can be es- 
timated using analytical approximations for the different 
asymptotic regimes ]22| ] . For the physical values used in 
Fig. [5] we find p* to be well within the rounded transition. 

Our simple model can also be applied to DNA melt- 
ing HQ, if the hydro gen bonds are considered to be 
generalized linkers. Existing models also show a sharp 
melting transition. The advantage of our approach is 
that some of the physical aspects of the chain interac- 
tions (namely, electrostatics and bending rigidity) are 
taken into account explicitly. On the other hand, the 
specific double helical structure of double stranded DNA 
is not included in the model. One straightforward gen- 
eralization of the model would be to include a sequence- 
dependent linking energy reflecting the different number 
of hydrogen bonds connecting adenine with thymine and 
cytosine with guanine. 

It is difficult to generalize our approach from two 
chains to many chains because many-body effects become 
important. Nevertheless, we expect our results to be re- 
flected in the behavior of many-chain solutions. We note 
that the "railway track" structure resembles a bundle 
where the chains are parallel and the linker density is 
high p4j . Similarly the dilute large- loop regime is remi- 
niscent of networks where the chains cross at large angles 
and the linker density is low. Such structures have been 
observed experimentally (4||| . Finally, in the limit of iso- 
lated chains the "railway track" is analogous to the com- 
pact torus formed in highly dilute DNA solutions jf],p|,^5| . 
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